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Abstract. We present a generalized version of the Path Interchange 
Zel'dovich Approximation (PIZA; Croft and Gaztahaga 1997; hereafter 
CG97), which may be used to reconstruct density and velocity fields from 
realistic galaxy redshift surveys. The original PIZA does not include the 
possibility of using galaxies in redshift space with a selection function. 
We map galaxy positions from redshift to real space and from the Local 
Group rest frame to the CMB rest frame. To take account of the selec- 
tion function, we minimize the mass-weighted action. We apply our new 
method to mock galaxy catalogues and find it offers an improvement in 
reconstructions over linear theory. Applying PIZA to the PSCz we can 
obtain real space density, initial density, and peculiar velocity fields, bulk 
velocities and the dipole. Comparison of our reconstructed dipole with 
the CMB dipole gives ~ 0.6. 



1. Introduction to PIZA 

CG97 noted that with the Zel'dovich Approximation, the least action solution 
is that which minimizes the mean square particle displacement, 



where the & are particle displacements. CG97 assume the galaxies have equal 
mass. The galaxies are assigned, at random, initial positions drawn from a 
homogeneous distribution. These are represented by PIZA particles. 

Applying PIZA is very simple. Two galaxies are picked at random and their 
particles are swapped. If S decreases, the swap is kept. This is repeated until S 
has been minimized. CG97 found that ~ 500 swaps per particle are needed. 

2. Problems with applying PIZA to redshift surveys 
2.1. PIZA in Redshift Space 

PIZA minimizes real space galaxy displacements, but we have redshift space 
displacements. In linear theory, the two are related by (Taylor and Valentine 
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where V^ 1 = 5fj — j^jSiSj is the inverse redshift space projection tensor. The 
displacement squared is then 
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2.2. Frames of Reference 

PIZA is carried out in the CMB frame, but the survey redshifts are in the Local 
Group frame. In the CMB Frame, 

i =t -tr + (1 + /?)2 • W 

To find a mock galaxy representing the Local Group is included in PIZA. 

2.3. The Selection Function 

In a flux limited survey, the number density of objects decreases with increasing 
redshift. We thus pick PIZA particles which obey the galaxy selection function, 
4>(r). To take account of galaxies below the flux limit, we assign each galaxy a 
mass of and each particle a mass of ^np, , where n is the particle-to-galaxy 
ratio. 



2.4. Sky Coverage 

There are regions on the sky, such as the Galactic Plane, that the survey does 
not cover. We pick PIZA particles with the same sky coverage as the survey and 
neglect PIZA trajectories that cross the plane. 



3. Two-Step PIZA 

We have developed a two-step PIZA method to include these generalizations. 
Two steps are needed because there are two possibly conflicting conditions to 
satisfy - the action must be minimized and the mass of each galaxy should be 
equal to the mass of its particles. 

In the first step, each galaxy is assigned n particles. The particle masses 
are used in equation (|l]), and interchanges are carried out as usual. 

The second step attempts to improve on step one, by getting M ga i axy ~ 
E m particies- A particle is taken from a galaxy with M galaxy < E ^particles (too 
many particles) and reassigned to a galaxy with M ga j axy > J2 m particles (too 
few). PIZA is then applied again, conserving the new particle to galaxy ratios. 
Generally this must be repeated 15 times. 



4. Tests on Simulations 

We have applied our new PIZA method to PSCz-like mock catalogues. We have 
tested the reconstruction of real space positions and radial peculiar velocities, 
and compared our reconstruction with linear theory. The mock catalogues and 
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Figure 1. Comparison of PIZA and Linear Theory Reconstructions. 
The left hand panel shows PIZA radial velocities compared with sim- 
ulation, the right hand shows linear theory compared with simulation. 
Velocities were smoothed with a Gaussian of radius 20 Mpc. The dot- 
ted line indicates the least-squares fit. 

linear theory reconstructions were provided by Enzo Branchini (Branchini et al 
1999). Our conclusions from these tests include the following: 

• The number of swaps needed is greater than that needed in the original 
PIZA - approximately 4 times as many. 

• Our reconstructions are generally better than linear theory (see Figure |l|) 

• The reconstructed Local Group velocity is not very sensitive to the as- 
sumed value of /3. 

5. PIZA applied to PSCz 

We have applied our new generalized PIZA method to the PSCz survey (Saun- 
ders et al 1995; also this volume), and present preliminary results here. Further 
results and analysis may be found in Valentine et al (1999). 

• The direction of the dipole was found to be ~ 20° away from the CMB 
dipole. Comparing the amplitude of these dipoles leads to a value of (3 ~ 0.6. 
Figure || shows the bulk flow and dipole. 

• We have compared our bulk flows with those of Mk III and Branchini et 
al, and find good agreement with both, implying a value of (3 ~ 0.7 

• The velocity field in the Supergalactic Plane is shown in Figure ||. 
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Figure 2. The PSCz bulk flow and dipole. The x-, y- and z- compo- 
nents and magnitude of the bulk flow and dipole are shown. It can be 
seen that there is no contribution to the dipole after 100/i _1 Mpc. 
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Figure 3. A slice through the PSCz velocity field containing the Su- 
pergalactic Plane. The slice is 20 fa _1 Mpc thick. The arrows show the 
galaxy displacements, and the arrow heads are the real space galaxy 
positions. 
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